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ABSTRACT 

Both the maximum size N m and the sea level muon size N y have 
been used separately to find the all-particle energy spectrum 
in the air shower domain. However the conversion required, 
whether from N m to E or from N y to E, has customarily been car- 
ried out by means of calculations based on an assumed cascade 
model. It is shown here that by combining present data on N m 
and N y spectra with data on 1) the energy spectrum of air show- 
er muons and 2) the average width of the electron profile, one 
can obtain empirical values of the N^ to E and N y to E conver- 
sion factors, and an empirical calorimetric all-particle spec- 
trum, in the energy range 2* 10 6 < E < 2* 10 9 GeV. 

1. Introduction. The great majority of shower particles are electrons, 
so it is natural that in the earliest air shower experiments the energy 
estimates were based on the number of electrons at the observation level. 
The first estimates (Auger 1939) were too low by about a factor of 4, be- 
cause 1) the correction for longitudinal development was too conserva- 
tive, and 2) the energy given to muons, neutrinos and low energy hadrons 
(E yV h) was ignored. (At the energies in question E yV h amounts to some 
35% of the whole.) The first difficulty stems from the average electron 
energy being comparatively low ('v E c , so that electrons are continually 
absorbed and regenerated. In order to estimate the energy deposited in 
the atmosphere one must use an integral signal such as the yield of at- 
mospheric Cerenkov or fluorescent light, or else face the problem of ac- 
curately evaluating a correction factor that may be as large as a factor 
20. By observing showers near maximum development (which generally means 
at a very high altitude) one can reduce the uncertainty in Ej;m by mini- 
mizing the correction factor. Following this approach, one finds the 
all-particle energy spectrum by combining measurements of the N m spectrum 
with estimates of the conversion factor E/N m . 

The alternative is to use shielded counters, which respond only to 
muons, and measure the N y spectrum, where N y (>E y ) is the 'muon size', the 
number of muons with enough energy to penetrate the shielding . This was 
done on a very large scale in the SUGAR experiment (Horton et al. 1983) 
and more recently in experiments at Chacaltaya, Tien Shan and Akeno 
(Kakimoto et al. 1981, Kirov et al. 1981, Hara et al . 1983). The diffi- 
culty with this method is that calculations relating N y to primary energy 
are relatively complex and model-dependent (see for example McComb et al . 
1977 and Hillas 1981) . Calculations of E/N m are less affected by these 
difficulties, but they also require estimating the energy given to muons. 

My purpose here is to show that by treating the experimental N m and 
N y spectra simultaneously , using also experimental data on 1) the energy 
spectrum of air shower muons and 2) the width of the electron profile, 
one can obtain conversion factors which are almost entirely empirically 
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based, and a new result on the all-particle energy spectrum which is al- 
most entirely model independent. 

2. Relation between maximum size and electronic energy. The energy dis- 
sipated by electrons is given by the track length integral 

e EM = (E c /x 0 )/N(x)dx , (1) 

where N(x) is the number of electrons at depth x g/cm 2 , E c is the criti- 
cal energy (=81 MeV in air according to Dovzhenko and Pomanskii 1964) , 
and x Q is the radiation length (= 37.1 g/cm 2 in air according to confer- 
ence paper HE4.4-4). Writing E EM = K (E c /x 0 ) aN m , where lf m is the height 
and a is the width (standard deviation) of the average shower profile, 
what can be said about the value of K? Using a Gaussian distribution for 
N (surely quite a crude approximation), K = /2i r = 2.51. Using a gamma 
distribution, N = N 0 £ q exp (-q£) where C == x/x m , which can be adjusted to fit 
very well (see conference paper HE4.4-5) , the value of K ranges from 2.35 
for q=6 (small showers) to 2.42 for q=12 (large showers). Thus it hardly 
varies at all, so adopting an average value for K, and substituting for 
(V, o) , X obtain Eem , (0/192) ^ , (2, 

where a is in g/cm 2 and E EM is in GeV, accurate to 1-2%.* The profile 
width has been measured in the Yakutsk and Utah experiments (Grigoriev et 
al. 1983, Baltrusaitis et al . 1985), but only for E ^ 10 9 GeV. The ener- 
gy dependence is expected on theoretical grounds to take the form a 2 = 

A + BD 10 logE, where D 10 is the elongation rate per decade, = 65 g/cm 2 
(Linsley and Watson 1981) and B is a characteristic length of order 60-70 
g/cm 2 (conference paper HE4.4-5) . Using the Yakutsk-Utah data to fix the 
value of A, one obtains a 2 = 1.1-10 4 + 4.2*10 3 logE (Linsley 1983),^ and 
finally, by substitution in (2) , 

E™ s 0.71 N™ 1 * 02 5 . (3) 


3. Relation between muon size and E^ vb . There is good agreement among 
independent measurements of the energy spectrum of air shower muons 
(Atrashkevich et al . 1983 and references therein). This spectrum is 
quite hard; almost half of the observed energy is given to particles with 
individual energies above 30 GeV. Over the range of shower energies 
where it has been studied (3»10 5 -10 8 GeV), the shape of this spectrum is 
invariant; hence the total energy of the observed muons is proportional 
to N u (>lGeV) , the number of muons (at sea level) with energy > 1 GeV, 
where the proportionality constant equals 10.0+.5 GeV. To obtain the 
energy given to neutrinos the observed muons are propagated backward to a 
production spectrum. In the air shower region it is found that E v ^ 

0.4 Ey f0 bs' where E v includes both and v e . This result checks with a 
forward propagation calculation by Hillas (1981) . Experiment based esti- 
mates of E^, the energy deposited by low energy hadrons, range from 
0.8 Ey ^ 0 bs (Greisen 1956) to 0.3 E^ obs . Adopting E h ^ 0.4 E y , obs as a 
conservative estimate, the total non-electronic contribution is obtained: 


E yvh “ < 18 GeV) -N (>lGeV) 


sea level 


An alternative form which may sometimes be more convenient is E EM = 
(Xhm/428) N m , where x^ is the full width at half maximum (Linsley 1981) 
* In the energy range of interest here, the simpler formula a = 130 + 
10. 2logE is equivalent. 
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Table 1. Muon size for a given in- 
tensity from various experiments • 
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Neither this result nor the one ex- 
pressed by (3) depends on any as- 
sumption about the primary composi- 
tion; they are properties of cosmic 
rays as they occur, in this energy 
range, at the solar system. As an 
experimental result, (4) applies to 
the energy range given above, 3 
10 8 GeV. Extrapolation up to 10 11 
GeV is justified unless there occurs 
some radical change affecting the 
production of very high energy muons 
and neutrinos. 


4. Calorimetric all-particle energy 
spectrum. Data on the Ny and N m 
spectra are summarized in Tables 1 
and 2, in inverse form (Ny and N m as 
functions of integral intensity) . 

The (inverse) all-particle energy 
spectrum is obtained by adding to- 
from (3) and Ey V ^ from 
tabulated values of Ny 
and for the same intensity, in 
the range icr 6 -10“ 12 /m 2 sr s where 
there are reliable data for both N 


gether E EM 
(4), using 
Nm 

range 


v 

and N m . The result, changed to dif- 
ferential form, is shown in Fig. 1. 

In case of Table 1 the various Ny 
values for a given intensity were 
averaged; in case of Table 2, the 
later results at Chacaltaya were 
hsed. There is good agreement with 
the Yakutsk energy spectrum in this 
range (Efimov and Sokurov 1983) , and 
with Haverah Park results (Cunning- 
ham et aim 1980). The present result 


supports a rather high location, nearly 10 7 GeV, for the transition re- 
gion where the change of slope (knee) occurs. 


5. Other results: conversion factors. As I showed previously, this deri- 

vation of the energy spectrum yields as by-products factors for convert- 

ing N m and separately to primary energy, and also yields the fraction 

of primary energy given to electrons, vs energy (Linsley 1983). These 

results are shown in Fig. 2. The low value found for E/N m , 1.3+.2, con- 

firms an important prediction by Hillas (1983) • An apparent conflict be- 

tween results from Chacaltaya and from lower elevations is resolved. The 

earlier values of N m (La Pointe et aim 1968) were somewhat too low, but 

energies were about correct because the conversion factor was somewhat 

too high. The later values of N m are more nearly correct, but the ener- 

gies were too high because the conversion factor was much too high. Al- 

though here the conversions, N m to E and Ny to E, have been treated sym- 

metrically, the energy dependence of the Ny-E 'conversion factor' makes 
this inconvenient in practice. Convenient formulae for representing the 
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N u data in Fig. 2 and Table 1 are: N y (>lGeV) s #1 = E (GeV) 0 * 8 3 s /6.8 , (5) 
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Figr. 1. All-particle Energy Spectrum . 



Fig m 2. Other Results . Open cir- 
cles E/N m , squares Egy/E ( r . h. 
scale) , filled circles logfE/N ^) • 


J[>N y (>lGeV)] s<1> - 3-10 4 N m 2 *\ (6) 

where J is in, m -2 sr -1 s -1 . 

Table 2, Maximum size for a given 
intensity from various experiments 
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